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SIMULATION MODEL OF A SINGLE-STAGE LITHIUM BROMIDE - WATER 

ABSORPTION COOLING UNIT 
by David Miao 
Lewis Research Center 

SUMMARY 

The performance and load capability of a given LiBr-HgO absorption chiller oper- 
ating with a hot -water heat source depends on six quantities: the inlet temperatures and 
flow rates of the hot -water source, the cooling -tower water, and the return chiller water. 
Based on this, a computer model for a single-stage absorption cooling machine has been 
developed which does not require data relative to the interior characteristics of the ma- 
chine (heat -transfer rates and surfaces). The model considers both heat -transfer and 
thermodynamic processes. It consists of two algorithms, one for the design, or refer- 
ence conditions, and the other for the off -design analysis. It is constructed from the 
steady -state equations but may also be used for the transient analysis of a cooling 
system. 

The program can be used in an independent mode or as a subroutine, as for example, 
with TRNSY’S, for the analysis of a cooling system. For a given size of machine the 
model can be used to predict off -design cooling -system performance, the only input re- 
quirements being a set of reference or rated conditions for the machine. 


INTRODUCTION 

The LiBr-HgO absorption liquid chiller has been used in the refrigeration and air- 
conditioning industry for some time. One of the primary reasons for using this type of 
machine is that steam or hot water, whichever is available, can be directly used as an 
energy source to power the machine. This characteristic is particularly attractive for 
solar -cooling applications. 

In a typical solar -cooling application, water heated through the passages of a bank 
of solar energy collectors is used to power an absorption machine to provide chilled 
water which in turn is used to air condition the building. 



Typically an absorption chiller is designed to handle the maximum expected load of 
the building. The design point thus represents a set of fixed operating conditions. How- 
ever, the actual load varies with building heat -transfer characteristics as well as local 
weather conditions. The design load may seldom be experienced. Since the chilled- 
water temperature is likely to increase with decreasing heat load (part -load operation), 
the chiller may be incapable of dehumidification. 

An approach is, prior to machine selection, to simulate various loading conditions 
through a computer model of the machine. The typical models available today are either 
empirical (ref. 1), or based upon a thermodynamic approach. The former generally 
represents a specific machine, and therefore its usefulness is limited; the latter is use- 
ful for providing a set of design conditions to the machine manufacturer to determine the 
size of an absorption machine. 

A thermodynamic approach can be used for simulating various operating conditions 
but such a model does not recognize limitations of the heat -transfer processes. A better 
approach is to take both heat -transfer and thermodynamic processes into consideration. 
Furthermore, if an existing machine is selected for a specific job, the heat -transfer sur- 
faces in the machine are fixed but often not known, and therefore it will be difficult to de- 
termine the capability of the machine over a range of operating conditions. The only 
known inputs are three sets of inlet flows and temperatures to the machine: namely, the 
flow rates and corresponding temperatures of the return chilled water, the cooling water, 
and the incoming hot water. The unknowns required to be established are the correspond- 
ing outlet temperatures of the three flow streams. A computer model for handling this 
type of problem is not generally available. The purpose of this report is to document a 
method for modeling the system. 


THERMODYNAMIC CYCLE 

The thermodynamic cycle of the absorption machine is well known (refs. 2 to 4). 
Figure 1 represents a typical arrangement of a single-stage machine. The machine bas- 
ically consists of five heat exchangers called a generator, a condenser, an evaporator, 
an absorber, and a solution heat exchanger. For a heat load imposed on the evaporator 
E, the LiBr-HgO strong solution is pumped through the solution heat exchanger X to the 
generator G. Heat energy is added in the generator G to drive out the refrigerant (in 
this case, water is the refrigerant). The remaining solution is called the weak solution. 
A portion of the weak solution is forced through the solution exchanger X and a pressure 
reducing valve VI back to the absorber A for the next cycle. To make a strong solu- 
tion in the absorber A, the refrigerant leaving the generator G must also be brought 
back to the absorber through a condenser -evaporator path. In the process, the refrig- 
erant is first condensed by removing its latent heat in the condenser C; then in passing 
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through an expansion valve V2, the pressure and the temperature of the refrigerant are 
reduced. The refrigerant is evaporated due to heat load addition in the evaporator E. 
The refrigerant vapor is then brought to the absorber A to be absorbed. When the heat 
of absorption is removed, the strong solution is restored and the new cycle begins. 

Most commercial machines are built on this basis. To simulate machine perform- 
ance, a thermodynamic cycle analysis is used to perform heat balance calculations in 
order to establish heat input and cooling requirement for a refrigeration load. Heat in- 
puts and outputs of the machine are marked in figure 1. The solid arrow lines indicate 
the direction of heat flow as well as fluid flow interior to the machine while the dashed 
lines indicate heat or fluid flowing into or out of the machine. Figure 1 is used to con- 
struct the thermodynamic portion of the machine model. 

To perform such calculations, thermodynamic properties of the water and the LiBr- 
HgO solution are also needed. Such information is readily found in reference 2 (the for- 
mulas used may be found in appendix A). 


HEAT -TRANSFER CONSIDERATIONS 

The thermodynamic analysis determines the cycle temperatures and the required 
heat flows for the five heat exchangers in the absorption machine as shown in figure 1. 

For a given refrigeration load, the heat exchangers must be designed to satisfy the afore- 
mentioned requirements. Once the heat exchangers are designed, the heat -transfer sur- 
faces are fixed and heat transfer is limited by the surfaces provided in the machine. 
Therefore, for all operating loads the performance of the machine is determined from 
the actual heat -transfer surface areas. 

In terms of heat -transfer processes on the LiBr-HgO solution side, the heat ex- 
changers may be classified into two types: The solution heat exchanger X (fig. 1) which 
deals strictly with sensible heat transfer is one type - and the other four: G, A, E, and 
C (fig. 1) which involve latent heat are another type. Heat exchangers G and A deal 
also with the heat of absorption, but since their heat -transfer coefficients are high and 
their temperature profiles are fairly constant, the heat -transfer analysis is treated in 
the same manner as those of exchangers E and C. The following equations (ref. 5) are 
used for these four exchangers: 


EFFN = 



( 1 ) 
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where 

EFFN temperature effectiveness of heat exchanger 
Tj inlet temperature of heating or cooling medium 
Tg outlet temperature of heating or cooling medium 

T temperature of LiBr-HgO solution of refrigerant (water) undergoing evaporation, 

absorption, or condensation process 

To relate the temperature field to the heat transfer, EFFN is rewritten as 

-UA/GC 

EFFN = 1 - e P (2) 

where 

U overall heat -transfer coefficient of heat exchanger 

A total heat -transfer surface 

G flow rate of heating or cooling medium 

C heat capacity of medium 

Jr 

Equations (1) and (2) are used to solve for the required outlet temperatures Tg 's 
of the four heat exchangers involving external fluid flows. Ideally, if all temperatures 
and flow rates are given at the design load, , equations (1) and (2) should resolve four 
UA’s for that machine. 

To simulate various heat loads other than the design, the corresponding UA's must 
be calculated from additional equations so that equations (1) and (2) can be used to obtain 
the various outlet temperatures Tg's. However, the information about the heat -transfer 
surface is usually not available and the UA terms are inseparable. Therefore, the 
next equations are derived on the UA term basis. 

Heat exchangers of this kind are typical shell -tube type. The cooling or heating 
medium is usually on the tube side, and the refrigerant (water) or LiBr-HgO solution is 
on the shell side. The heat -transfer process is governed by the mechanism of the fluid 
flow on both sides and the tube wall thermal resistance. By definition U is written as 

I-J.+I + R+F (3a) 

u N, » 

where 

h w tube -side coefficient due to forced convection 

h shell -side coefficient 
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Rj. tube wall resistance 
F sum of fouling factors on both sides 

The term in formula (3a) is a function of the tube wall thickness and the mate- 
rial of construction. Typically, its magnitude is very small because of low pressure 
operation and the use of high conductivity copper based tube material. 

The design or selected fouling factor F (ref. 6) is also rather small. The true foul- 
ing factor varies with water conditioning and plant operation and cannot be established 
without test data. Both R^ and F may be considered constant throughout machine 
operation. 

The h term, due to latent heat transfer, is very high for a good cost effective heat 
exchanger design. The h value for boiling water or steam condensation may be on the 
order of two to six times the forced convection coefficient h^ (ref. 7). Therefore, it is 
not a strong factor on the overall heat -transfer coefficient U, which may be conveniently 
written as 


where R is the sum of the resistances (1/h) + R^. + F. 

Equation (3b) implies that U can be found if h^, is known. 

To find on the tube side, the following forced convection formula for turbulent 
flow (ref. 7) is used: 

(4 a 

where 

D inside diameter of tube 

K thermal conductivity 

H viscosity 

A flow area 

V 

Equation (4a) indicates that the change of h^ is sensitive to the changes of the flow 
rate G (eight -tenth power function) but less dependent on the heat transport properties. 
Furthermore, the fluid temperature variations for an absorption machine are rather 
small, especially in a solar application; thus these temperature dependent properties 
remain practically constant. Therefore equation (4a) may be rewritten as 
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Since proportionality can be established, may be written as follows in terms of 
a reference condition with the subscript 0: 

hw = (g~) hw0 <4b> 

ts) 

By combining equations (3b), (3c), and (4b) and solving for U, we obtain 


U = 


/_G_\°- 8 ( 1 + h wO R p \ 
V G o/ V 1+h w R / 


u f 


(3d) 


As long as the term h^R is not substantially different from h^R^, 4116 * actor 
(1 + h w0 R 0 )/(l + l^R) is approaching unity. If this term is assumed to be one, the ex- 
pected error in U is 5 to 10 percent. Under the worst conditions, the error may be as 
high as 20 percent. Therefore, equation (3d) may be reduced to 


or 


U = 



U 


0 


UA = 8 (UA) q (3e) 

Equation (3e) implies that, if a reference condition is known, the UA term at other 
conditions can be found given the right flow proportions. To find a reference UA, equa- 
tions (1) and (2) must be used and flow rates are referred to the reference condition. Us- 
ing actual measured values in the aforementioned formulas instead of the machine design 
values for the reference point is desirable wherever possible. 

The second type of heat exchanger in the absorption machine is a liquid to liquid ex- 
changer (exchanger X in fig. 1). This exchanger is placed in the absorption circuit to 
improve cycle efficiency. It is also typical of a shell -tube type with a true counterflow 
arrangement for better heat recovery. The strong solution (rich with water refrigerant) 
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is pumped through the tubes and the weak solution flows across the tube bundles, with 
flow deflecting baffles. As was pointed out previously, the heat -transfer rate is a strong 
function of the flow rate. The strong solution flow rate is greater than that of the weak 
one. To achieve a high heat -transfer coefficient on the tube side, it is natural for the 
heat exchanger designer to place the strong solution in the tubes. In addition, the better 
heat transport properties of the strong solution (more water content) aid in achieving a 
high coefficient. The lower shell -side coefficient of the weak solution can be improved 
by using spaced baffles. 

Equation (4a) is used to calculate the tube-side coefficient. Equation (4b) is also ap- 
plicable if the heat transport properties remain practically constant. 

As indicated previously, equation (4a) or (4b) is applicable for turbulent flow. For 
a true counterflow type of heat exchanger, or single-tube pass arrangement, the velocity 
in the tubes may be reduced under some part load operation. It is possible the flow pat- 
tern may shift into the laminar region. Then equation (4a) or (4b) would not be appli- 
cable, and the formula for laminar flow (ref. 7) would have to be used. 

Since this report is concerned with the simulation of a previously designed machine 
without knowing the interior arrangement of the heat -transfer surface areas, the laminar 
formula, even if it is available, is probably not useful for model construction. How- 
ever, it is reasonable to assume that the turbulent flow formula is used for calculating 
the tube -side heat-transfer coefficient. In these machines, the heat exchanger with 
longer tube lengths (thus small flow area and high velocity in the tube) is commonly seen 
in commercial machines. 

The formula for the shell-side coefficient (ref. 7) may be written as follows because 
the heat transport properties remain practically constant: 


WV.0.33 

K 


D G 
e w 


0.6 




cross; 



or 


where 


h or G 6 
gW W 


(5) 


D e equivalent diameter 

h coefficient of weak solution flow rate 

gw 

A„„ noo flow passage area measured along shell inside diameter 
cross 

Unlike the tube -side formula, equation (5) is not restricted by the turbulent flow. 
The shell -side coefficient can be increased by means of closer baffle spacings. 
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Therefore, it is reasonable to assume that the weak solution with less flow rate is on the 
shell side. 

The relation between the overall heat transfer and the individual coefficients is the 
same as that of equation (3a). In this case the controlling resistance is on the tube side 
because of the single tube pass arrangement. The magnitude may be on the order of the 
shell -side coefficient. Since heat -transfer coefficients on both sides are poor, the mag- 
nitude of (l/h w ) + (1/h) in equation (3a) is much larger than that of and F (perhaps 
10 times larger); therefore, R^. and F are neglected and equation (3a) may be re- 
written as 



where 

x refers to solution exchanger 
gs refers to strong solution 
gw refers to weak solution 

For a referenced condition, equation (6a) becomes 


(6a) 


1 


U. 


xO 



(6b) 


Once again for a given machine, where the interior construction of the machine is 
not known, equation (6b) cannot be solved without making assumptions. If hg g g and 
kgwO are assume d equal, equation (6b) becomes 


h gs0 " h gw0 ~ 2U x0 


( 7 ) 


By combining equations (7) and (4a) or (5), h and h can be obtained for other sim- 

gb gw 

ulated conditions; specifically 






( 8 ) 


( 9 ) 
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Then substituting equations (8) and (9) into equation (6a) and rearranging the terms yield 


U x = < 2U xO> 



( 10 ) 


Since the heat -transfer surface area is fixed, equation (10) may be written as 



Equation (11a) again shows that the overall heat -transfer rate at any other condition 
can be established through a known reference condition (design or test). Equations (8) 
and (9) can also be extended to include the property corrections if better accuracy is de- 
sired. The heat transport properties except thermal conductivity may be found in ref- 
erence 3. For thermal conductivity values for various LiBr-HgO solutions, a fraction of 
water conductivity proportional to water concentration are suggested. In general these 
effects on heat -transfer coefficients are small and will not be taken into consideration at 
this time. 

The aforementioned equations were derived on the assumption that hg S Q = hgwO’ *" he 
assumption appears valid because (1) the fluid properties on both shell and tube sides are 
similar and (2) the flow rates are not substantially different within the operating range of 
the solution concentration. However, if hg S g is substantially different from h^Q, 
equation (11a) may be generalized as 


(UA) X = (Fl)(UA) x0 



(lib) 
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where, for example, 


h gsO h gwO 
h gsO « h gwO 
h gsO = X ‘ 5h w0 


FI = 2 and F2 = 1 


FI = 1 and F2 = 0 


FI = 2. 5 and F2 = 2/3 


Equation (lib) may be useful to experimentally determine the actual values of FI 
and F2 for use in the program for a given machine. 

Equation (11a) or (lib) can be solved if (UA) x q is known or may be found from a 
given set of the design temperatures. The effectiveness is given in terms of the tem- 
peratures (refs. 1 and 4) as 


T - T 

EFFNX = _i (12) 

T - T 

g a 


where 

T temperature of generator 

S 

Tg outlet temper ature of weak solution 

T temperature of absorber 
a 

In general the exchanger is designed with the effectiveness EFFNXq = 0. 7 to 0. 8. 
If Tg in equation (12) for the design load is not known, the relation between EFFNXq 
and T 50 may be established by heat balance (ref. 1). 

When EFFNXq is found together with flow rates G g Q and G w q and the solution 
heat capacities C g Q and C w q, then (UA) x q is calculated from the following equation: 


-NTU x [(l-(C mln /C max )] 

EFFNX = — 

1 /^min \ ^x [^"^min^max^] 


(13) 


where 

NTU (UA) /C . 
x v V min 

C . G C 
min w w 

C G C 
max ss 
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C w heat capacity of weak solution 

C_ heat capacity of strong solution 
s 

The subscript 0 used previously has been deliberately omitted in equations (12) and 
(13) for the purpose of generalization. Then C min = (G W C W ) and C max = (G g C g ) be- 
cause G„. < G„ and C < C„ for LiBr-H 0 0 absorption machine. The (UA)„ n is solved 

w S W S ct aU 

implicitly in equation (13). 


OTHER CONSIDERATIONS 

The equations derived in the previous section together with the thermodynamic equa- 
tions discussed in the section THERMODYNAMIC CYCLE are the working formulas for 
the five heat exchangers to be used in the construction of the simulation model. In addi- 
tion to these formulas, heat losses, pump capacity, operating range of the solution con- 
centrations, and operating temperature limits should be included. Unfortunately machine 
construction does vary with the design approach of different manufacturers, and the con- 
struction information is usually not available. It is difficult to generalize all the limit- 
ations to be accommodated by the model. Nevertheless, some of the important consider- 
ations that should be taken into account follow. 


Heat Losses 

The heat losses vary with the specific design and the ambient environment in which 
the machine is installed. Heat may leak out of or into the machine, and between the 
partition shells separating the heat exchangers in the machine. The result is that addi- 
tional heat supply is required to accommodate these losses. To account for these losses, 
a simplistic approach is to add a fixed percentage to the heat supply. A few percent may 
be sufficient for the type machine considered herein. The thermodynamic equations 
(appendix A) may be modified as follows: 


«G - < G w H 5 - G s H l + GrH,)(FGQ) 


Qg - g r(h 7 - h 8 ) 


Qp 

1 + 1— — Id - f qc) 


(Qq + Q' 


E; 


Q a - (G w H 5 - G s Hj + g r h 10 ) 


1-1 — -Hi-Fqo) 


\Qq + Q 


E / 


(14) 

(15) 

(16) 
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where 



multiplication factor 
1 (no heat loss considered) 

1. 02 (equivalent 2 percent loss) 


Solution Pump Capacity 

Normally the pump capacity is chosen to meet the design load. For part load oper- 
ation, the required flow rate may or may not exceed the maximum capacity. For a par- 
ticular load demand, if heat source temperature is low and/or the cooling water temper- 
ature is high, the machine, based on the thermodynamic cycle analysis, tends to demand 
more solution flow. Since the flow control is not known and varies somewhat with differ- 
ent machines, it is assumed that the solution flow rate cannot exceed the capacity of the 
design point. 


Concentration of the LiBr-HgO Solution 

For an absorption process to exist in operation, there are limits on the solution con- 
centrations. If the concentration is too rich, crystallization will occur. If the concen- 
tration is too lean, no absorption process will occur. Reference 2 suggested that the 
concentrations should be kept within 0. 5 to 0. 65 range. For this model a range from 
0.4 to 0. 68 has been used. 


Temperatures and Temperature Differences 

The temperature limitations, like the solution concentrations, are set for the oper- 
able absorption process. Usually these are the outlet temperatures of the external fluids 
in heat exchanger G, C, A, and E (fig. 1). The limits of these temperatures have been 
placed in the program (see appendix B). 

In addition to the temperature limits, the temperature difference across the heat ex- 
changer surfaces are also limited by the heat -transfer processes. In general the tem- 
perature differences between the two heat exchange mediums at outlet condition will be 
used for setting the limits (see appendix B). 

When the aforementioned limits and the concentration limits are properly set, the 
solution heat exchanger temperature as well as the pressure limits may be neglected. 
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MODELING ALGORITHM 


With the necessary equations and the limiting conditions established, the next step is 
to formulate an algorithm for computer operation. The desired solution for a given set 
of inputs is the one that achieves the lowest possible outlet temperature of the chilled 
water. The heat balance is not only required to satisfy the thermodynamic analysis but 
also simultaneously satisfy the heat exchanger equations. 

The model consists of two different algorithms. One part is used to solve for the 
reference or design conditions. Another part is used to solve for the off -design condition 
based on the established reference condition. The second part is simply to perform an 
internal heat balance to establish the corresponding outlet temperatures of the three flow 
streams, namely hot water GH, cooling water GC, and chilled water GE. The calcula- 
tion sequence for this part is first outlined as follows: 

(1) Input GH, GC, GE, TH1, TA1, TE1 and an off -design tonnage, (see fig. 1). 

(2) Calculate flow rate per ton for flow GH, GC, and GE. 

(3) Calculate effectiveness (eq. (2)) for exchanger G, C, A, and E. 

(4) Calculate TE2, TE, TH2, and TG. 

(5) Calculate TC2 with an assumed COP. 

(6) Assume TA. 

(7) Calculate TC. 

(8) If TA or TC exceed limits, change tonnage. 

(9) Calculate TG, TC, TA, and TE with newly assumed tonnage. 

(10) Calculate solution concentration. 

(11) If XI or X4 exceeds limits, change tonnage. 

(12) Calculate enthalpies H8 and H10 of refrigerant at outlets of condenser C 

and evaporator E, respectively. 

(13) Calculate refrigerant flow GR and solution flows GS and GW, respectively. 

(14) Calculate effectiveness EFFNX of solution exchanger. 

(15) Calculate two outlet temperatures T3 and T5 of solution exchanger. 

(16) Calculate refrigerant enthalpy H7 at outlet of generator G, weak solution 

enthalpy H5 at outlet of solution exchanger X, and strong solution enthalpy 

H, at outlet of absorber A. 

(17) Calculate generator heat QG, condenser heat QC, and absorber heat QA. 

(18) Calculate COP. 

(19) If TA is not agreeable with assumed value, adjust TA to suit. 

(20) If COP is not agreeable with assumed value, adjust COP to suit. 

(21) Check temperature difference limits. 

(22) Check pumping rate limits. 

(23) Check concentration limits. 

(24) Force tonnage to maximum. 



(25) Check chilled water outlet temperature TE2 at set point. 

(26) Calculate pressure PE and PC. 

To establish the reference conditions, several of the aforementioned indicated steps 
are repeated. The algorithm used depends upon the information available. 

If all the design or reference temperatures are given but the flow rates are not, 
steps (10) to (18) and step (26) are repeated. The flow rates and all reference (UA)’s 
are the calculated outputs. The effectiveness of the solution heat exchanger can be cal- 
culated from the known temperatures (eq. (12) as an input to the program). 

If all three external flow rates are known instead of their outlet temperatures, 
steps (1) to (18) and step (26) are repeated. In this case the corresponding outlet tem- 
peratures are determined. 

If the outlet temperature of the solution heat exchanger T5 or its effectiveness is 
not given, an assumed effectiveness must be used as an input until a rated reference ton- 
nage is found. 


PROGRAM DESCRIPTION 

The computer program was written in FORTRAN IV language. It can be used as a 
subroutine to simulate the absorption machine performance in a cooling system. Al- 
though the equations derived are steady -state type, no restriction is imposed for use in 
the transient analysis of a cooling system. 

When used as a subroutine, the program may have to be modified to accept a set of 
the design or the test conditions. The flow rates and the inlet and outlet temperatures of 
the three external fluid streams are system connected to run the simulation. If addi- 
tional outputs such as heat loads, COP, and operating pressures are required, they may 
be system linked or printed out for analysis. 

When used as an independent program, the first case is treated as the reference 
case. The program calculates additionally needed reference values and stores these 
values automatically in the program. Starting with the next case, the user inputs as 
many off -design cases as are desired. NAMELIST input is used in the program. 

All tolerances for the limitation conditions discussed previously have been prestored 
in the program but can be changed as desired. The units system used to perform the cal- 
culation is metric but provision to use British units for inputs and outputs is included. 
Changing either the units or the tolerances shall be discussed in the next section and 
appendix B. 
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OPERATION OF THE PROGRAM 


Use as a Subroutine 

If the program is used as a subroutine, the reference data and program controls 
must be inserted as data statements or their equivalent by the user. The required data 
are UAGO, UACO, UAEO, UAAO, UAXO, GSO, GWO, and TONO. The controls are FQG, 
METRIC, KLBHR, and JWRITE (see appendixes B and C). 

The input variables are currently placed in an array called XIN. These variables 
(listed in order), are GHT, GCT, GET, TH1, GA1, TE1, and TONX (see appendixes B 
and C). 

The output variables are arranged in an array OUT. These variables are GHT, 
GCT, GET, TH2, TC2, TE2, and TON. If additional outputs such as COP, PC, and PE 
are required, the user may place these variables in the additional locations of array 
OUT (see appendixes B and C). 


Use as the Main Program 

If the program is used independently, the reference data must be calculated from 
this program based upon the available design or experimental informations. The input 
variables in this case will be TH2, TA2, or TC2, TE2, TH1, TA1, TE1, TG, TC, TA, 
TE, TONO, TONX, KLBHR, METRIC, and JREF (see appendix B and fig. 1). TONX is 
the initial guess of the actual load. The data are entered via a NAMELIST read and are 
for reference case. The NAMELIST name is REF. The first tabulated output will be the 
results of the design conditions and the table is identified with a case marked 0. 

To run other cases with fixed heat -transfer surfaces (the same machine), additional 
cases are placed in the run stream with a NAMELIST name of VAR. As many cases as 
desired can be run. The input for these cases are GH, CC, GE, TH1, TA1, TE1, and 
TONX (see appendixes B and C). The outputs are tabulated as before, and the case is 
identified with a case number greater than 0. 

The convergence is controlled by KTA, KCOP, KTONI, and KTON2. If the number 
of the iterations is excessive, the output may be incorrect. The user must examine the 
results to decide whether he should increase the number of iterations, or discontinue his 
rim because of exceeding machine operating constraints. 

The tolerance controls for the temperatures and concentrations are currently pre- 
stored in the program (see appendixes B and C). The values may be changed to suit the 
user's purpose. 
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SAMPLE CALCULATIONS 


Two sample computer printouts are included to demonstrate the use of the program 
in appendix C. 

Sample 1 shows that, for a given set of the design conditions, the program not only 
finds the correct design load but generates the results for the off -design loads as well. 

The absorption machine used in the sample calculations is a TRANE model C1H 
(ref. 6). This model was designed for a nominal rated tonnage at 174 tons. The print- 
out table (case 0) shows that the calculated tonnages agree with the design load. The out- 
put of this case is then stored in the program as the reference data of the machine to be 
used for the off -design rims. 

A total of 130 off -design cases (the off -design loads and operating conditions in 
table 2C1H of ref. 6), have been run with the program. Most of the calculated tonnages 
agree with the data in reference 6 within 2 percent and generally are slightly greater 
than the table values (two typical cases are shown in appendix C). In some of the cases, 
however, the calculated values are high by 9 percent. These cases usually are asso- 
ciated with the extremely high or low outlet temperature of the chilled water. All cases 
were run on the assumption that the nominal design flow rates were chosen to establish 
the rated table values. If these flow rates are not nominal but varied within the design 
range, the program calculated tonnages can be brought to agreement with those tables 
indicated. 

Sample 2 was intended to show that, with minor changes, the program can be used 
as a subroutine in a system program. In this case the system program is TRNSYS 
(ref. 1). Sample 2 is a solar assisted building cooling system modeled with TRNSYS 
program (see appendix C). 


CONCLUDING REMARKS 

A computer model of a LiBr-HgO single-stage absorption machine has been devel- 
oped. By utilizing a given set of design data but without knowing the interior character- * 

istics of the machine, the off -design performance can be simulated or evaluated. Al- 
though the model is not validated experimentally, it can be a useful tool for analyzing the V 

capability of a given machine, or for studying the machine performance in a cooling 
system. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 16, 1978, 

776-22. 
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APPENDIX A 


SYMBOL LIST, THERMODYNAMIC FORMULAS, 
AND EQUATIONS FROM REFERENCE 1 

Strong concentration (XI > 0. 5): 


Xi _ (49 .04 + 1. 125 TA - TE) 
(134. 65 + 0.47 TA) 


kg LiBr 
kg solution 


Weak concentration (X4 < 0. 65): 


X4 = (49.04 + 1. 125 TG - TC) 
(134.65 + 0.47 TG) 


kg LiBr 
kg solution 


Enthalpy of condenser outlet: 


H8 = (TC - 25) kcal/kg 


Enthalpy of evaporator outlet: 


H10 = (572. 8 + 0. 417 TE) kcalAg 


Refrigerant flow: 


GR = 


QE kg 


(H10 - H8) hr 


Strong solution flow: 


GS = GR 


X4 kg 


(X4 - XI) hr 


Weak solution flow: 


GN = GR *1. 

(X4 - XI) hr 
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Heat capacity of strong solution: 


CXI = 1. 01 - 1. 23(X1) + 0. 48(X1) 2 kcal/(kg)(°C) 
Heat capacity of weak solution: 

CX4 = 1. 01 - 1. 23(X4) + 0.48(X4) 2 kcal/(kg)(°C) 
Outlet temperature of weak solution 

T5 = TG - (EFFNX)(TG - TA) °C 
Outlet temperature of strong solution: 


T3 = T A + (EFFNX) / — \ \ (TG - TA) °C 

\ x4 / \CX1) 


Enthalpy of absorber outlet: 

HI - [42. 81 - 425. 92 (Xl) + 404. 67(X1) 2 ] + [l.Ol - 1.23(X1) + 0. 48(X1) 2 ](TA) 


Enthalpy of weak solution at heat exchanger outlet: 

H5 = [42.81 - 425. 92 (X4) + 404. 67(X4) 2 ] + [l.Ol - 1. 23(X4) + 0. 48(X4) 2 ](T5) 


Enthalpy of refrigerant at generator outlet: 

H7 = (572. 8 + 0. 46 TG - 0. 043 TC) kcalAg 

Condenser heat load: 

QC = (GR)(H7 - H8) kcal/hr 

Generator heat load: 

QG = (GW)(H5) + (GR)(H7) - (GS)(H1) kcal/hr 


kcal/kg 


kcal/kg 


t 

1 
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Absorber heat load: 


QA = (GW)(H5) + (GR)(H10) - (GS)(H1) kcal/hr 
Coefficient of performance: 




cop = 

QG 


♦ Evaporator heat load: 


QE = 3024. 0 kcal/hr 


Evaporator pressure: 


PE = antilogy 7. 8553 - 


1555 11.2414x10^ 


TE + 273. 15 ( T e + 273.15) 2 


mm Hg 


Condenser pressure: 


PC = antilog 10 7. 8553 - 


1555 11.2414X10^ 


TC + 273. 15 ( TC + 273 15) 2 


mm Hg 


* 

4 
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APPENDIX B 


SYMBOL LIST FOR HE AT -TRANSFER CALCULATIONS IN COMPUTER PROGRAM 
Flow rates, gal/min, lb /hr, kg/hr 


GH 

(Hot water supply) 


GC 

(Cooling water supply) 


GA 

(Cooling Water supply) 


GE 

(Returning chilled water) 

i 

GR 

(Refrigerant - water) 


GW 

(Weak solution) 


GS 

(Strong solution) 


Temperatures, °F, °C 
TH1, TH2 

(Inlet and outlet conditions of GH) 


TCI, TC12, or TA2 

(Inlet and outlet conditions of GA) 


TC12, TC2 

(Inlet and outlet conditions of GC) 


TE1, TE2 

(Inlet and outlet conditions of GE) 


TG 

(Generator) 


TC 

(Condenser) 


TA 

(Absorber) 


TE 

(Evaporator) 


Heat -transfer rates, Btu/(hr)(°F), cal/(hr)(°C) 


UAG 

(Generator) 


UAC 

(Condenser) 

* 

UAA 

(Absorber) 

V 

UAE 

(Evaporator) 


UAX 

(Heat exchanger) 


First digit 

Overall heat -transfer coefficient 


Second digit 

Overall heat -transfer surface 


Third digit 

Component symbol 
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Number of heat -transfer units 


NTUG 

(Generator) 

NTUC 

(Condenser) 

NTUA 

(Absorber) 

NTUE 

(Evaporator) 

NT TJX 

(Heat exchanger) 

Heat -transfer 

effectiveness 

EFFNG 

(Generator) 

EFFNC 

(Condenser) 

EFFNA 

(Absorber) 

EFFNE 

(Evaporator) 

EFFNX 

(Heat exchanger) 


(A digit 0 following aforementioned symbols signifies a reference or a design condition 
being used. A digit T following aforementioned symbols and symbols in appendix A 
signifies total quantities. ) 


TONO 

(Reference refrigerant tonnage) 

TON 

(Tonnage calculated) 

TONX 

(Tonnage variable) 

COPX 

(COP variable) 

TAX 

(TA variable) 

GSC1 

(Product of strong solution flow and heat capacity) 

GWC4 

(Product of weak solution flow and heat capacity) 

CRATIO 

=GWC4/GSC4 

EXPX 

(Exponential function for heat exchanger) 

Controls and 

limits 

METRIC 

(Input to be metric units > 0) 

KLBHR 

(Input to be lb/hr > 0) 

JWRITE 

(Write output > 0) 

KTA 

(TA converging cycle = 50) 

KCOP 

(COP converging cycle = 50) 
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KTON1 and KTON2 
ACONST = 1. 0° C 
BCONST = 1. 296° C 
CCONST = 1.425° C 
DCONST = 1. 919° C 
TELO - 2.22° C 
TE2SET = 4.43° C 
COPHI = 0. 93 
COPLO = 0. 60 
FQG = 1. 0 
EFFNX = 0.71428 


(TONX converging cycle = 100) 
Limits of (TE2 - TE) 

Limits of (TA - TC12) 

Limits of (TC - TC2) 

Limits of (TH1 - TG) 

(Lowest temperature limits of TE) 
(Lowest temperature limits of TE2) 
(Highest limits of COP) 

(Lowest limits of COP) 

(No heat loss added) 

(Initialization of EFFNX) 



APPENDIX C 

SAMPLES 1 AND 2 WITH PROGRAM LISTINGS 


l* 
2 * 
3* 
4* 
5# 
6* 
7* 
8 * 
9* 
10 * 
11 * 
12 * 
13 * 
14* 
1 5* 
16* 
17 * 
19* 

1 9* 
20 * 
21 * 
22 * 
23* 
24* 
25* 
26* 
27* 

2 8 * 
29* 
30* 
31* 
32* 
33* 
34* 
35* 
36* 
37* 
38* 
39* 
40* 
4 1* 
4 2* 
4 3* 
4 4* 
4 5* 
4 6* 
4 7* 

4 8* 
49* 
50* 
51* 
52* 
S3* 

5 4* 
5 5* 
56* 
57* 
58* 
$9* 
60* 
61* 
62* 
63* 


Sample 1: LiBr-^O Single-Stage Absorption Machine Used as a Main Program 

C USE THIS TO EVALUATE OUTPUT OF AN ABS0PPTI0N MACHINE WITH FIXED --UA — 

C ALL WATEP SPECIFIC HEAT £ DENSITY ASSUMED TO BE —1.0— EXCEPT HOT WATER 

DIMENSION XIN(10),PAR(13),XNTU(6),EFFN(6) 

DIMEN SI ON X<6)*Y(6),GIN<3> 

DIMENSION T0NGVN<16P),TONCAL<160) 

C 

C--METR IC-0 , BRITISH UNITS USED. JWRITE=1 WRITE ALLtJWRITE-C NO WRITE 

C KLBHR-OfGPM FOR FLOW INPUT. KLBHR=l f LBS/HR INPUT 

DATA METRIC/O/ , KLBHR/Q/ , JWRITF/1/ 

DATA PDKC/.4536/ 

DATA TFTCl/32./,TFTC2/1.8/ 

DATA CALBTU/3. 96831/ 

C 

C 

C CONST 1 £ CONST 4 ARE CONCENTRATION LIMITS 
DATA CONST1/0.4/.CONST4/C. 68/ 

C A-B-C-D-CONST ARE LIMITS FOR EVAP., ABSOPP., COND.,£ GENERATOR 
DATA ACONST / 1 . /,BC0NST/1. 296/.CC0NST/ 1.423/ tDCONST/1 .919/ 

C 

DATA TELO/2 .22/ , T E 2S E T/4. 4 3/ 

C COP LIMITS HEAT LOSS FACTOR 

DATA COPHI /D.9 3/,C0PL0/n. 6f / , fqg/ 1 .0/ 

C EFFNXrD. 71426 FOR T5=135 F E F F NX - ( TG-T 5 ) / ( T G -T A ) 

DATA EFFNX/O. 71428/ 

C 

DIMENSION X I N R S V ( ID ) 

DIMENSION A J T E 2 ( 4 ) 

DATA A JTE2 < 1 > , A JTE ?< 2 ) , A JTE2 ( 3 > , AJTE2 < 4 > / 3HTH1 , 3HTE1 , 

1 3HTri? , 3HTE2 / 

DIMENSION A JREF < 2 > 

DATA AJPEF ( 1 ) ♦ A JREF ( 2 )/ 3HTC2 , 3 H T A 2 / 

DATA JTE2 / G / 

DTLIM-C.25 

DTLIMrC.l 

NAMELIST /REF/ TH2 , T A 2 , TC 2 f TE? t TH 1 , TA 1 , TE1 » TG , T C , T A , TE f TONE,, 

1 KLBHR, METRIC, JREF, CO NS T1,C0NST4,AC0NST,BC0NST,CCCNST,DC0NST, 

2 C0PHT,C0PL0,FQG,EFFNX,TEL0,TE2SET 

NAMELIST /VAR/ GH,GC,GE,THl,TAl,TEl,T0NY,JTE2,TE2,TH2tLVAP 
C 

C JPEFrO.RUN NO REF. -JREF-1,TC2-XIN(2) • — JREF-2,TA2 OR TC12 = XIN(2) 

3 E D READ( 5 ,REF ) 

L V A R -O 
XIN i 1 > = TH2 

I F ( JREF .EQ . 1 ) XIN(?)-TC2 

IF<JREF.E0.2) XIN(2)-TA2 

X IN ( 3 ) =TE2 

XINC4 ) - T H 1 

X IN I 5 ) =TA 1 

T C 1 - T A 1 

XIN(6)~TE1 

X I N ( 7 ) = TG 

X IN < 8 ) ~ T C 

XI N ( 9 ) =TA 

X I N < 1 D ) - T E 

DG 3 □? I- 1 , 1C 

X IN R S V ( I ) - X IN < I ) 

3 P 2 CONTINUE 

T C2RSV-TC? 

T A2RSV-T A2 
KERROP-C 
C 
C 
C 
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64* IF< METRIC. GT.L ) TFTCl=n.D 

65* IFC METRIC-GT,0) TFTC2=1.0 

66* IF(METRIC.GT.C) PDKG=1.0 

67* IF( METRIC. GT*0) CALBTU-1.0 

68* IF(HETRTC«GT.C ) BPH=l.r 

69* C 

70* BPH=5DO.G 

71* IF< KLBHR.GT.D ) BPHrl,0 

72* C 

73* IF< JREF-GT.O) GO TO 19 

74* 7 CONTINUE 

75* C UA VALUES ARE PER TON BASIS 

76* TONO-TON 

77* G SD-G S*P DKG 

78* G WO-G W*PDK G 

79* GHO=GHT 

80* GEO-GET 

81* GCOrGCT 

82* P AR < 1 )3UAG*PDKG 

83* PAR(2)-UAC*PDKG 

84* PAR(3I-UAE*PDKG 

85* P A R ( 4 ) -UA A*PDK G 

86* PAR(S>=UAX*PDKG 

87* U AXO-P AP ( 5 ) 

88* C 

89* C 

90* INPUT - 1 

91* C 

92* 99D CONTINUE 

93* JREF-0 

94* READ( 5, VAR t END-999 ) 

95* XIN(1 )-GH 

96* X I N ( 2 > - GC 

97* X I N ( 3 ) - GE 

98* X IN < 4 ) “ TH 1 

99* XI N ( 5 ) — T A 1 

100* T C 1 - T A 1 

101* X IN ( 6 ) “ TE 1 

102* X I N ( 7 ) - TON X 

103* GESAV-XIN 13) 

104* TE2 SA V- TE 2 

105* TH2SAV-TH? 

13 6* K ERRO P-Q 

107* C TON IS AN ASSUMED VALUE TO ST A D T 

108* T ON X- X I N ( 7 ) 

109* TONREF=TONO*D .5 

110* T ON MI N-TONO*0« 1 

111* T ONMAX-TONC*! -2 

11 2* IFC XIN(7> .LE* TONREF ) TONXrTONREF 

113* C 

114* 19 CONTINUE 

115* C 

116* K GC -G 

117* 9 CONTINUE 

118* C 

119* C FLOW PATES ARE PER TON BASIS 

120* C PEP TON BASIS 

121* GE- 33 23-9573 

122* W 1 - 0 * 

123* R 5 - D • 

124* H7=C. 

125* H8=0. 

126* H 1 0 - 0 • 

127* TONrQ.O 

128* K TO N 1 -0 

129* K T 0 N2 -0 

130* KT0NX--1 » 
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1 31* 
1 32* 
133* 
134 * 
135* 
136 * 
137* 
138* 
139* 
14 3* 
141* 
142* 
143* 
144* 
145* 
146* 
147* 
148* 
149* 
150* 
151* 
152* 
153* 
154* 
155* 
156* 
157* 
158* 
159* 
1 6 G* 
161* 
1 62* 
163* 
164* 
165* 
166* 
167* 
168* 
169* 
170* 
171* 
172* 
173* 
174* 
175* 
176* 
177* 
178* 
179* 
1 R0* 
161* 
182* 
163* 
184* 
185* 
186* 
187* 
188* 
189* 
1 9 C* 
191* 
192* 
193* 
194* 
195* 
196* 
197* 


X4=D.67 
X 1=D .41 
I WRI TF = 0 
11 CONTINUE 
C C0PX-.722 

IF ( CO P X # L F • COPLO#OR»COPX*GE*COPHI ) C0PX~D*722 

K C 0 P -0 
16 CONTINUE 

GH = XIN< 1> /T0NX*PDKG*6P H*D • 9 75 
GC-XIN (2) /T0NX*PDK6*GPH 
GF-XINC3) /TONX*PDKG*BPH 
THlr(XlN(4)-TFTCl)/TFTC2 
TC1=(XIN(5>-TFTC1>/TFTC2 
TE 1 - ( XIN(6 )-*TFTCl)/TFTC2 
T E 2 S V Mr | TF2SAV-TFTC 1 )/TFTC2 
TH2SVN=(TH2SAV-TFTC1)/TFTC2 
C 

IF( JRFF.EO.O) GO TO 20 
TH2=(XIN(1)-TFTC1)/TFTC2 
T A 2 - ( X I N ( 2J-TFTC1 )/TFTC2 
IF( JREF.EQ.2) TC12-TA2 
IF< JRFF.EO.l) T C 2- T A 2 
TE2-{XIN( 3)-TFTCl)/TFTC2 
TG=< X I N ( 7 ) -T FT C 1 J/TFTC2 
TC-(XIN(8)-TFTC1)/TFTC2 
T A — < XIN(9)-TFTC1)/TFTC2 
T E - ( XINC 10)“TFTC1)/TFTC2 
GO TO 21 
20 CONTI NUF 

C ASSUME TUBE SIDE WATER FILM COEF. IS CONTROLLING 

GH1 = C X I N ( 1 ) /G HO ) ** 0 • 8 
3 C 1 - ( XI N ( 2 ) /G CO )**3 .8 
G El - ( XIN( 3)/GFC)**0*8 
C 

C TOTAL BASIS IN METRIC UNITS 

GIN ( 1 ) -Gh / G H 1 *T 0 NX 
GIN(2)’GC/GC1*T0NX 
G I N I 3 )-GE/GEl*TONX 
C 

DO 10 I “ 1 f 4 

X N T U ( I ) -PAR ( I ) /GIN ( I >* 1 • 0 * TO NO 
IF ( I . E0.4 ) XNTU(I)-PAR(I)/GIN(2)*T0N0 
IF ( XN T U ( I ) • GF . IT * ) GO TO 8 
EFFMI)=1.C-EXP(-XNTUII )) 

GO TO 10 
6 CONTINUE 

EFFN( I>rD.999 
1C CONTINUE 
C 

c 

TE2-TE1-QE/GE 
TH2=THl-( OE/COPX)/GH 
C 

C FOR CHECKING TRANE TABLE FIGURES ONLY TE 2 + TH2 KNOWN INSTEAD 

IFIJTF2.LF. C ) GO TO 603 
TE2-TE1 
TH2-TH1 

C T E2 “ T F2 S A V 

C T H2“ TH2 S A V 

TE2-TF2SVM 
TH2=TH2SV* 

TE1=TE2+QF /GE 
THl ” T H 2 + { OE/COPX ) / GH 
6CC CONTINUE 
C 
C 

T E- TE 1 - ( TE1-TE2 )/FFFN< 3 ) 
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198* 
199* 
230* 
201 * 
202 * 
203* 
2D4* 
205* 
206* 
207* 
2D 8* 
209* 
210 * 
211 * 
212 * 
213* 
214* 
215* 
216* 
217* 
216* 
219* 
220 * 
221 * 
222 * 
223* 
224* 
225* 
226* 
227* 
228* 
229* 

2 3D* 
231* 
232* 
233* 
234* 
235* 
236* 

2 3 7* 
238* 
239* 
240* 
24 1* 
242* 
24 3* 
244* 
245* 
246* 
247* 
248* 
249* 
250* 
251* 
252* 
253* 
254* 
255* 
256* 
257* 
258* 
259* 
260* 
261* 
262* 
263* 
264* 


TGrTHl-(TNl-TH2)/EFFNt 1 ) 

TC2-TC1+C 1.0+ l.O/COPX )*QE/GC 
C ASSUMED A VALUE FOP TA 
TA2-(TC1+TC2)*0.5 
TA3TC1-(TC1-TA2)/EFFN C 4 ) 

K T A - 0 

15 CONTINUE 

TC-TC2/EFFN(2 )-(l .0/EFFN(2 > -1 . 0 > * C TC 1 +EFFN ( 4 > * < T A -TC 1 >) 

IFC TC.LF.TA ) GO TO 4G 
IF(TE.GE.TA ) GO TO 41 
IF(TC.GF.TG) GO TO 41 
C 

21 CONTINUE 

XI = ( 4 9.04 +1. 12 5*TA-TE) /( 1 34 . 65+D . 47*T A ) 

X4= <49.04+1. 1 2 5* T G -TC ) /(I 34 . 6 5 +D . 47* TG ) 

I F ( X 1 .LT.C0NST1 ) GO TO 45 
IF(X4.LE.X1) GO TO 43 
C 

H6-TC-25.0 
H10-572.6+0.417*TE 
GR=OE /( H10-H6 ) 

G$rGR*X4/ CX4-X1 ) 

GW=GS*< XI /X4 ) 

C 

CXlr 1.0 1-1.2 3*X 1+0. 48 *X 1**2 
CX4-1 *0 1*1 • 23*X4+0.48*X4**2 
G SC 1 “ GS *CX 1 
5 WC4rGW*CX4 
CRATI0-GWC4/GSC 1 
I F f JREF.GT.O) GO TO 22 
C 

C ASSUMING ORIGINAL FILM COEF. ECU AL ON BOTH SIDES. --GW ON SHELL SIDE 
Flr2.D 
F2rl.O 

C Fl=2 # F2=l*HGS = HGU.~Fl = l,F2 = 0,HGS>>HGW.~Fl = 2*5 f F2 = 2/3,HGS = 1.5HGW. FOR UAX 
RGS - ( ( GSC/GS) * ( TONG/TONX) )**J.8 
RGU- < <GWD/GWI*< TO ND/TON X i ) ** j . 6 
UAX-F1*UAX0*( 1 .0/ (RGS+F2*RGW) ) 

C 

IF! GWC4 .GT.GSC1 ) GO TO 17 
XNTUX-UAX/GWC4*(T0N3/T0NX ) 

GO TO 18 
C 

17 CONTINUE 

XNTUX-UAX/GSC 1* (TONQ/TONX ) 

CRATI0=GSC1/GWC4 

18 CONTI N UF 
C 

r F( ABS( 1. n-CRATIO ) .LT .0.01 ) GO TO 13 
IF< <XNTUX*( 1.-CRATI0) J.GE. 10. > GO TO 12 
EXPX-EXP( -XNTUX*( l.O-CRATIO) ) 

EFFNX-( l.H-EXPX)/( 1 .D-CRATIO*FXPX ) 

GO TO 14 
C 

12 CONTINUE 
EFFNXZO .999 
GO TO 14 

13 CONTINUE 

E FFNXrxN TUX/11. O+XNTUX! 

14 CONTINUE 
22 CONTINUE 

C 

C 

T5rTG-EFFNX* ( TG-TA ) 

T3 = TA+( EFFNX*CRATIO*( TG-TA ) ) 

C 

Hl:(42.81-425.92*X1+4D4.67*X1**2)+CX1*TA 
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265* 
266* 
267* 
268* 
269* 
2 7 D* 
271* 
272* 
273* 
274* 
275* 
276* 
277* 
278* 
279* 
2 8D * 
281* 
282* 
283* 
284* 
285* 
286* 
287* 
288* 
289* 
290* 
291* 
2 92* 
293* 

29 4* 

2 9 5* 
296* 
297* 
298* 
299* 
300* 

30 1* 
302* 
333* 
304* 
305* 
306* 
337* 
300* 
309* 

3 1 D * 
311* 
312* 
313* 
314* 
315* 
3 16* 
317* 

31 8* 
319* 
320* 
321* 
322* 
323* 
324* 
325* 
326* 
327* 
328* 
329* 
330* 
331* 


H5-(42.81-4 25.92*X4+4 04 .67*X4**2> + CX4*T5 
H7-572.8+n.46*TG-0.043*TC 
C 

C F 0 G- 1 • 0 » NO HEAT L OS S . — F GG > 1 .0 , L <2 . 0 , X HEAT LOSS 

IFCFQG.LE. 1.0) FQG-1 .0 
QG-< 6W*H5-GS*H1+6R*H7 )*FQG 
QCrGR*< H7-H8) * < 1 . G + Q G / < QG + QE ) * < 1 . 0 -F& G ) ) 

GA=< GW*HS-GS*H1+GR*H1D >*t 1 •0 + Q6/< QG+QE )*< i .Q-FQG) ) 
C OP -Q E / QG 
C 

IF( JREF.GT.O) GO TO 60 
C 

T C12-TC1 +QA/GC 
T AX-TC1-ITC1-TC12 )/EFFN(4 ) 

C 

40 CONTINUE 

I F < TC.LE.TA ) TAX-TC 

IF(ABS(TAX-TA).LT. 0.00001) 60 TO 41 

IF* KTA.EQ.5D) GO TO 41 
T A- ( T A X + T A ) * 0 • 5 
K T A T A + 1 
GO TO 15 

41 CONTINUE 
C 

I F < A 6 S ( COPX-COP ).LT .0*00001) GO TO 42 
IFCKC0P.EG.50) GO T042 
COPXr { COPX + COP ) *D . 5 
K COP-K COP + 1 
GO TO 16 

42 CONTINUE 
C 

X ( 1 ) -TE2-TE 
X(2)-TA-TC12 
X { 3 > - T C- T C 2 
X (4 ) z TH2 -T G 
Y ( 1 )~ ACONST 
Y (2 ) - B C ONS T 
Y ( 3 ) -CCONST 
Y < 4 ) "DCONST 


DO 47 1-1,4 

IF( X< I) .LT . YC I ) . AND .KTONX .EO. 1 ) GO TO 60 
I F ( X ( I ) * L T • Y I I ) ) GO TO 45 
47 CONTINUE 
C 

IF(X1. GT. C0NST1.AND.X4.LT. C0NST4. AND. X4.GT. XI) GO TO 46 
C 

45 CONTINUE 

I F ( K T ON 2-100) 49,43,43 

49 CONTINUE 

T ON MI N - T ON X 

T0NX-(T0NX+T0NMAX)*0.5 
I TONX-IFIX (TONX) 

TCNXrTONX + OTLTM 
KT0N2-KT0N2+1 
GO TO 11 
C 

43 CONTINUE 

I F ( KT0N1-100) 44,65,65 

44 CONTINUE 

T ONHAX-TONX 

T ONX-C TONX+TONMIN )*0.5 
I TONXrIFIX( TONX ) 

IF< TON. LE. 0.0. AND. TONX. LE. (TONMIN+l.G) ) GO TO 60 
IF< TONX. LE. TON ) GO TO 5D 

TONX-FLOAT ( I T 0 N X ) - 1 .0 
TONX-TONX - DTLIM 
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K T 0 N 1 - K TON 1 + 1 
GO TO 11 


332* 

333* 

334* 

335* 

336* 

3 3 7 * 

338* 

339* 

340* 

341* 

342* 

343* 

344* 

345* 

346* 

347* 

348* 

349* 

350* 

351* 

352* 

353* 

354* 

355* 

356* 

357* 

358* 

359* 

363* 

361* 

362* 

3 6 3* 

364* 

365* 

366* 

367* 

368* 

369* 

370* 

371* 

372* 

373* 

374* 

375* 

376* 

377* 

378* 

379* 

380* 

381* 

382* 

383* 

384* 

385* 

386* 

387* 

388* 

389* 

3 9D * 
391* 
392* 
393* 
394* 
395* 
396* 
397* 
398* 


C 

46 CONTINUE 

IFCKTONX.EQ.l.OR.TONX.LE.D .0) GO TO 60 
IF(TON.GE.TONX) GO TO 5D 

C CHECK MAX. STRONG SOLUTION PUMP RATE 

GSTG-G$P*TONC 

GSPUMP-GS*TONX 

IMGSPUMP.GT.GSTO) GO TO 43 
C 

T ON- TON X 

IMTE2.LT.TE2SET.0R.TE.LE.TELO) GO TO 48 
GO TO 49 
48 CONTINUE 

TON-TON*( TE1ME2SCT ) / (TE1-TE2) 

C 

50 CONTINUE 

IF(KTONX.EQ.I.OR.TON.LE.Q.O) GO TO 60 
T ON X -TON 
K TON X r 1 
GO TO 11 
C 

60 CONTINUE 

65 IF (JREF.GT.O.OR.KTONX.EQ.l .OR.KGC.GE .30) GO TO 66 
K6C-KGC+1 

X I N ( 3 )-XIN(3)*(TONX-1.0)/T ONX 
GO TO 9 

66 CONTINUE 
C 

C 

C ALL SPEC HEAT = 1. 

I F( JREF.EO. 0) GO TO 35 
T ON-TONO 

YY = ( 1 .-EFFNX*CRATIO )/ ( l.-EFFNX ) 

XNTUX-ALOG(YY )/( l.-CRATIO ) 

GH= Q G / ( TH1-TH2) 

GEr QE/(TE1-TE2> 

IF (JREF.EO. 1) GC-(QA+QC>/ (TC7-TC1 ) 

I F C JREF.EO. 2) GC = 0 A / ( T C 1 2 -TC 1 ) 

IF(JREF.EO.l) TC 12-TC1+QA/GC 
IF( JRFF.LQ.2 ) TC2-TC 12+OC/GC 
EFFNU ) = ( TH1-TH2) /(TH1-T6 ) 

EFFN( 2)=(TC12“TC2 )/ (TC12-TC) 

EFFN( 3)-(TEl“TE2)/(TEl-TE ) 
EFFN(4)r(TCl-TC12)/(TCl-TA) 

DO 31 K - 1 * 4 

XNTU(K)“(“*1. )*ALOG( l.-EFFN (K) ) 

31 CONTINUE 
C 

35 CONTINUE 

U AG=XNTU( 1 )*GH 
UAC-XNTU( 2 ) *GC 
UAE-XNTUt 3 ) *GE 
U AA-XNTU( 4 )*GC 
U AXrxNTUX*GWC4 
C 

ANALOG (10.0) 

9-1555. D/ (TE+273.15) 

C-11.2414E4/(TE+273.15)**2 
PE = EXP( A*(7.8 553-E-C) ) 

P-1555. D/(TC+273. 15) 

C-l 1.241454/ ( TC+273. 15)**2 
PC = EXP( A* (7.8 55 3*9-0 ) 

C 

OG=QG*C ALPTU 
0 C-Q C*C AL BTU 
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399* 


QE=QE*CALBTU 

400* 


QArQA*CALBTU 

431* 


H 1-H1*C ALBTU 

402* 


H 5-H5*C AL BTU 

403* 


H 7 “H7 * C AL BTU 

404* 


H 8 - H8 *C A L R T U 

435* 


H1D~H1D*CALBTU 

406* 

C 


407* 


T3=T3*TFTC2*TFTC1 

408* 


T5-T5*TFTC2+TFTC1 

439* 

C 


410* 


TH1=TH1*TFTC2+TFTC1 

411* 


TH2=TH2*TFTC2+TFTC1 

412* 


TC1=TC1*TFTC2+TFTC1 

413* 


TC12=TC12*TFTC2+TFTC1 

414* 


T A2-TC12 

4 1 ; 


TC2=TC2*TFTC2+TFTC1 

416* 


TE1=TE1*TFTC2+TFTC1 

417* 


TE2=TF2*TFTC2+TFTC1 

4 18* 


TE = TE *TFTC2+TFTC1 

419* 


TA =TA *TFTC2+ TFTC 1 

423* 


TC - T C *T FTC 2 +T FTC 1 

421* 


TG =TG *TFTC2+TFTC1 

422* 

c 


423* 


U AG=U AG/PDKG 

424* 


U AC- U AC /P DK G 

425* 


UAE=UAE/PDKG 

426* 


U AArU A A/PDKG 

427* 


U AXrUAX/PHKG 

428* 


G H- GH /P DK 6 

429* 


GC-GC/PDKG 

430* 


G A- GC 

431* 


GE-GE/POKG 

4 32* 


GR-GR/PTKG 

433* 


GW-GW/PDKG 

434* 


GS-GS/PDKG 

435* 


GSC1-GSC1 /PDKG 

436* 


GWC4-GWC4/PDKG 

437* 

c 


439* 


3GT- QG*T0N 

439* 


qct=qc*tom 

443* 


Q ET -T ON 

44 1* 


3 AT ~q a*tqn 

442* 


U A G T “ U A G*T ON 

443* 


U AC T -UA C* TON 

444* 


UAET-UAE*TON 

445* 


UAAT-UA A*TON 

446* 


GHT-GH*T0N/BPH/0 *975 

447* 


GCT-GC*TGM/BPH 

448* 


G A T - G C T 

449* 


GET=GF*TON/BPH 

453* 


DT12E=TE 1-TE2 

451* 


DT12A-TA2-TC1 

452* 


DT12C-TC2-T A2 

45 3* 


DT126-TH1-TH2 

454* 


DTE2=TE2-TE 

455* 


0TA2-TA-TA2 

456* 


DTC2-TC-T C2 

457* 


D TG2 - T H? - T G 

458* 

c 


459* 


IF<KTA.GE.5C.OR.KCOP.GE.5G> GO 

460* 


IF< KT0N2.GE . 1 00 . OR . K TON 1 • G E • 1 0 0 ) 

461* 


IFC JWRITE.EO.Ol GO TO 59 

462* 


GO TO 430 

463* 


58 CONTINUE 

464* 


K EPROPr l 

465* 

c 
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466* 

467* 

468* 

469* 

47D* 

471* 

472* 

473* 

474* 

475* 

476* 

477* 

478* 

479* 

483* 

48 1* 
482* 
483* 
484* 
485* 
486* 
487* 
488* 
489* 
490* 
491* 
492* 

49 3* 
494* 
495* 
496* 
497* 
498* 
499* 
500* 
531* 
502* 
503* 
504* 
5D5* 
536* 
507* 
508* 
509* 
510* 
511* 
512* 
5 13* 
514* 
515* 
516* 
51 7* 
516* 
519* 
520* 
521* 
522* 
523* 
524* 
525* 
526* 
527* 
528* 
529* 
530* 
531* 
532* 


430 I F ( METRI C. GT. 0) WRITE(6 t 420) 

420 FORMAT! lHl,20X, ’OUTPUT IS IN SI UNITS DFGREES C, KG/HP, CAL * / ) 

IF(METRIC.EQ.O) WRITE (6*421) 

421 FORMAT ( 1H 1 , 20 X , f GUTPUT IS IN (J.S. CUSTOMARY UNITS - DEGREES Ft GPM 

1, BTU • / ) 

IF(KERROR.EG. 1) WRITE(6,431> 

431 FORMAT! // / 20X , *********** ITERATIONS FAILED TO CONVERGE ******** 
1 ** • ) 

I F ( METRI C.EO* 0 .AND. INPUT .EO.O) WRITE(6,93> A JR E F ( JP E F ) 

IF! METRIC. GT.C .AND. INPUT. EQ. 3) WRITF(6,193) A JR E F ! JP EF ) 

IF! INPUT *E0.0) WRITE (6,94 ) C X IN til , I = 1 , 1 0 > , TONO , JRE F , INPUT 
IF! INPUT. EO.O) GO TO 57 
J OU T “ 1 

IF ( JTE2.6T.0 ) JOU T - 3 

IF(METRIC.GT.C) WRITE(6 t 195) AJTE2 (JOUT) ,AJTE2 (JOUT+1 ) 

IFC METRIC. EO.O) WRITE(6,95> AJTE2(J0UT),AJTE2(J0UT+1) 

IF! JTE2.E0.0) 

lWRITE(6 f 98) KTA t KC0P,KT0N2 ,KTONl f C XIN Cl ) , I rl ,6 ) , TONX ? XIN(7 ) ,INPlJT 
IF! JTE2.GT .0) WRITE <6, 98 > KT A , K COP , KTON? , K TON 1 , XIN ( 1 > , X IN ( 2 > , 

1 XIN(3),TH2SAV,XIN(5),TE2SAV,T0NX,XIN(7>,INPUT 
57 CONTINUE 


C 


WRITE (6 f 403) 
URI TE t 6 , 402 ) 
WRITE<6,4P7) 
WRITE! 6 , 4 P2 ) 
WRITE ( 6,406 ) 
WRITE! 6,405 ) 
XT ,QET 

WRITE(6,404) 
XT ,QAT 

WRITE <6,401 ) 
XCT 


X1,X4,CX1,CX4,GR,GS,GW,GSC1,GWC4,EXPX,CRATI0,UAX 
TA,T5,T3,TG,H1,H5,H7,H8,H10,XNTUX,EFFNX,C0P 
GE , T E 1 , TE 2 , TE , XNTU ( 3 ) , E FFN ( 3 ) ,PE,UAE,QE,GET,UAE 

G A , T C 1 , T A 2 ,TA , X NTU ( 4 ) ,EFFN( 4), PE, UAA,QA, GAT ,UAA 

GC,TA2,TC2,TC ,XNTU(2) ,EFFN(2),PC*UAC,0C,GCT,UACT t 0 


WRITE(6,4D0 ) GH , T H 1 
XT ,QGT 

4 C 0 FORMAT! IX, 4 HG ,7F10.3,5E1D.3//) 

401 FORMAT! IX, 4 HC ,7F10.3,5E10.3/> 

402 FORMAT! 1X,4HX v 12F10«3/> 

403 FORMAT! 5X , 12DH XI X4 

X GS GW GSC 1 GWC4 

X ) 

404 FORMAT! IX, 4 HA ,7F10.3,5E10.3/> 

40 5 FORM AT! IX, 4 HE , 7F10.3,5E 10.3/ ) 

406 FORMAT { 5X , 12DH G 

X EF FN P 

X ) 

407 FORMAT! 5X, 120H TA 

X H 5 H 7 

X ) 


,TH2,T6 ,XNTU(1),EFFN(1),PC*UAG,QG,6HT,UAG 


CXI 

EXPX 


CX4 

CRATIO 


GR 


U A X 



T 1 

T 2 

T 

NTU 

UA 

0 

GT 

UAT 

OT 


T5 

T 3 

TG 

HI 

H6 

H10 

NTUX 

effnx 

COP 


WRITE(6,96) 

WRITE<6,97> DT12E,DTE2,TE2,DT12A,0TA2,TA2,DT12C,DTC2,TC2,DT12G,DTG 
X 2 , TH 2 

9 3 FORMAT! 1 OX, •TH2-F’ , 5X,A3,•-F’,5X,’TE2-F*,5X,*TH^-F*,5X,’TC1-F , , 

1 5X, f TE1-F* ,5X, •TG-F* f 6X f ’TC-F»,6X,’TA-F’,6X,’TE-F’,4X,’T0N-PEr’, 

2 5X , • JREF * , 5X , ’NO . • ) 

193 FORMAT! 10X, *TH2^C’,5X,A 3, ♦ -C ’ ,SX, • TE 2-C ’ ,5X, ’ T H 1- C * ,5X , ' T C 1- C * , 

1 5X , * TE1 -C * ,5X ,*TG-C’ ,6X, *TC-C * , 6X , # TA-C ’ , 6 X , # TE -C ’ , 4 X , * TON-REF * , 

2 5X , * JREF’ ,5X, ’NO. * ) 

94 FORMAT! 5X , 1 1FI0.3 *110, 3X f I 3/ > 

95 FORMAT! 1UX, *KTA KCOP KT0N2 K T ON 1 GHT-GPM CCT-6P 

1 M GFT-GPM* ,5X,A3,’-F TC 1 -F • , 5 X , A 3 , • -F TON-CAL TON-START 

2 NO.* ) 

195 FORMAT! 10X, ’KTA KCOP KT0N2 KT0N1 GHT-KPH GCT-KP 

1H GET-KPH’, 5X, A3, ’ -C T C 1 - C ’ , 5 X , A 3 , ’ - C TON-CAL TON-START 

2 NO • * ) 

96 FORMAT! EX, 12CH DT12E DTE2 TE2 DT 1 2 A D T A2 
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533* 



X T A 2 DT12C DTC2 TC2 

534* 



X ) 

535* 


97 

fORMAT(lX,4HTEMP f 12FlG.3///> 

536* 


96 

F0RMAT(5X f 4(I6 f ?X) f 8F10.3f3X f I3/> 

537* 

C 



538* 


59 

CONTINUE 

539* 



I F < JREF.GT. 0) GO TO 7 

540* 

C 



541* 


6 1 

CONTINUE 

542* 

C 



543* 



T ONGVN ( INPUT ) -X IN ( 7 ) 

544* 



TGNICAL < INPUT ) = TON 

545* 



I NPUT = 1 ♦ INPUT 

546* 



G H-X I NCI ) 

547* 



G C - X I N ( 2 ) 

548* 



G E - GE S A V 

549* 



TH1-XTN(4) 

550* 



T A 1 -X I N ( 5 ) 

551* 



TE1-XINC6 ) 

552* 



T0NX=XIN(7) 

553* 



TE2=TF2SAV 

554* 



TH2=TH2S AV 

555* 



IF(LVAR.EO.O) GO TO 990 

556* 



DO 301 1=1,10 

557* 



XlNRSVf I)=XINCX) 

558* 


301 

CONTI NUE 

559* 



TC2=TC2RS V 

560* 



T A2=T A2RSV 

561* 



GO TO 3no 

562* 


999 

CONTINUE 

563* 



WRI TE (6 , 601 ) 

564* 

c 


WRITE (6,602) (TONGVN(M) t TONCAL (M ) ,M- 

565* 



WRITE<6,602) (M,TONGVN ( M ) * TON C AL ( M ) f 

566* 


601 

FORMAT ( lXf 11DHTNE FOLLOWING APE KNOWN 

567* 



1 MODEL C1H/DS-ARS1 TH2,TE2,TC1 ARE 

568* 

c 

602 

F0RMAT(lX,6(4X f ?F8.2)/ /(lX f 6<4X,2F6. 

569* 


6C2 

F0RMAT(lX,6(I4 f 2F6.2)/ /(lX t 6CI4,2F8. 

570* 

c 



571* 



STOP 

572* 



END 


DT12G 


DTG2 


T H 2 


TRANE 

) 




OUTPUT 

IS IN U.S 

. customary 

UNITS - 

DEGREES F, 

, GPH, BTU 







TM2-F 
233 .330 

T A2-F 
95.003 

TE2-F 
44 .000 

TH1-F 

270.000 

TC1-F 

85.000 

TE1-F 
54 .000 

TG-F 
210. GOO 

TC-F 
112 . 00C 

T A-F 
105.000 

TE-F 

40.000 

TON-REF 
174. OHO 

JREF 

2 

X 

XI 

.587 

X4 

. 643 

CXI 

.453 

CX4 

.420 

GO 

12.C07 

GS 

145.916 

GW 

133.909 

GS C 1 
66.161 

GUC4 

56.197 

EXPX 

.000 

CRATIO 

.849 

U AX 

119 .237 

X — 

TA 

IC5 .D00 

T5 

135.001 

T 3 

168.705 

TG 

210.000 - 

HI 

1 95.964 

H5 

-158.932 

H7 

2445.978 

He 

77.162 

HIP 

2280.403 

NTUX 

2.122 

EFFNX 

.714 

COP 

.721 

E 

G 

1199.983 

T 1 

54.000 

T 2 

44 .000 

T 

40.000 

NTU 
1 .253 

EFFN 

.714 

P 

6.234 

UA 

, 150*04 

Q 

.120*05 

GT 

.418*03 

DAT 

. 262*06 

OT 

.174*03 

A--- 

1573 .675 

85 .003 

95.0PD 

105.003 

.693 

.500 

6.234 

. 109+04 

.157+05 

+548*03 

.190*06 

.274*07 

O — 

1573.675 

95.000 

103.1*8 

112. ODO 

• 658 

.482 

69.925 

. 1G4 +04 

.129*05 

.548*03 

. 180+06 

.224+07 

G 

415.964 

270. COD 

230. 0P0 

210 .003 

1 .099 

.667 

69.925 

.457+3 3 

.166*05 

• 1 48*03 

.795*05 

.290+07 

TEMP 

DT12E 

13.330 

DTE2 
4. 003 

TE2 

44.000 

DT 12A 
ID .033 

DT A 2 
10.033 

T A ? 

95 .DOD 

DT12C 

8.198 

DTC 2 
8 .832 

TC2 

1D3 . 198 

0T12G 
40. OOC 

DTG2 

20.000 

T H2 

230. OOC 




OUTPUT 

IS IN U.S 

. CUSTOMA 

RY UNITS - 

0EGRFES F, 

, GPH, BTU 







X TA 
18 

K C CP 
IS 

KT0N2 

□ 

KTON 1 
0 

GHT-GPH 

150.033 

GCT-GPM 

553.000 

GET-GPM 
41 7.6 n 0 

TH 1-F 
273 • 47D 

TC1-F 
80 .ODD 

TE1-F 

5C.4Z8 

TON-CAL 
161 .902 

TON-START 
181 .OCO 

X— 

X 1 

.588 

X4 

.651 

CXI 

.453 

CX4 

.413 

GP 

1 1 .978 

GS 

123.773 

GU 

111.794 

GSC 1 
56.04 5 

GWC4 

46.147 

FXP X 
.671 

CRATIO 

.823 

U A X 

104.353 

X 

T A 

130 .341 

T5 

129 . 187 

T 3 

166.356 

TG 

2T9 . 364 

HI 

-200.533 

H5 

-161.478 

H7 

2445.731 

H 8 

67.912 

H 1 C 

2276.534 

NT UX 
2.261 

EFFNX 

.735 

C Ur 
.734 

E 

G 

1147. 87U 

T 1 

50. 428 

T? 

39.974 

35.792 

NTU 

1.253 

EFFN 

.714 

P 

5.276 

UA 

.144*04 

0 

. 1 2 n +C 5 

GT 

.4 18+03 

OAT 

•2b2+r6 

OT 

.182+03 

A 

1520. 046 

80. ODD 

90.157 

ino. 34i 

.692 

.499 

5.276 

. 1 D 5 + 3 4 

. 1 54 + 0 5 

.553+r3 

. 191 +[ 6 

• 2 8 1 +U 7 

c 

1523 .346 

93 . 157 

98.656 

107.804 

. o5 7 

.482 

61.929 

.999+03 

.129+05 

.553+03 

. 1 82 + 06 

• 25 5 +C 7 

G 

4C2.DD2 

273.470 

229 .779 

2 P9 . 364 

1.096 

.666 

61.929 

.44 1 +03 

.164+05 

.15r+P3 

. 30 2+1' 5 

.29g+I 7 

TEMP 

P T 1 2E 
10.454 

r TE 2 
4 . 182 

TE? 

39,974 

DT 12A 
1C . 157 

DT A 2 
10.184 

T A 2 

90.157 

DT 12 C 
8.499 

DT C 2 
9.149 

TC2 

98.656 

D T 1 2G 
40.691 

DTG2 
2P.4 15 

TH 2 

229. 779 




OUTPUT 

IS IN U.S 

. CUSTOMA 

RY UNITS - 

DEGREES F 

, GPM, BTU 







KTA 

K C OP 

KT0N2 

K T ON 1 

GHT-GPM 

gct-gpm 

GET-GPM 

TH1-F 

TC 1-F 

TE 1 -F 

TON-CAL 

TON-ST ART 


18 

1 7 

2 

2 

150.000 

553 .000 

417. 6 n D 

265.096 

BC.oao 

5C.G81 

175.864 

172. OOC 


X 1 

X4 

CX 1 

CX 4 

G D 

GS 

GU 

GS C 1 

G WC 4 

EXPX 

CRATIO 

U AX 

X 

.585 

.645 

.455 

.4 16 

1 1.966 

128.347 

116.38 1 

58. 338 

48.422 

.684 

. 8 30 

IDS .254 


TA 

T5 

T 3 

TG 

HI 

H 5 

H7 

H8 

H 10 

NTUX 

EFFNX 

CCP 

X 

99.670 

128 .275 

1 64.246 

206.075 

-201.447 

-163 . 724 

2442.486 

65. 827 

2276 . 662 

2.236 

.731 

.734 


G 

T 1 

T 2 

T 

NTU 

EFFN 

P 

UA 

0 

GT 

U A T 

OT 

E 

1187.279 

50 . 081 

39.974 

35.931 

1.253 

.714 

5.3P6 

. 149+34 

.120*05 

•418*n3 

• 262+P6 

• 17fc+r3 

A 

1572.235 

80.000 

89.822 

99 .67D 

.6*2 

.499 

5*306 

. 109+04 

.154+05 

.553*03 

.191*06 

.272*07 

c — 

1572.235 

89. 822 

98.027 

106.859 

.657 

.482 

60.239 

. 1 D 3 +D 4 

.129+rs 

« S53+P3 

. 18?*rfc 

• 22 7+r 7 

G 

415.804 

265. G9& 

225.7*3 

206-075 

1.0*6 

. 666 

60.239 

.456+03 

. 163+05 

• 1 5P+03 

.802+05 

.287*07 


DT12E 

DTE2 

T E2 

DT 1 2 A 

D T A 2 

T A 2 

DT12C 

DTC2 

TC2 

DT 1 2 G 

DTG2 

TH 2 

TEMP 

13 .137 

4.043 

39.974 

9.822 

9.848 

89.822 

8.205 

8.832 

93.027 

39.303 

19.718 

225.793 


Case 

U 


Case 

1 


Case 

z 
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I 


Sample 2: LiBr-HgO Single-Stage Absorption Machine Used as a Subroutine in TRNSYS 


1* SUB & OUT I NT T Y PE 1 7 ( T I M F t XIN t OUT * T ,0107 »PAR ) 

2* COMMON /PR2/ TIME? *TF INAL*DELT 

2* C USE THIS TO F V AL U A T E OUTPUT OF AN ABSORPTION MACHINE WITH FIXED — UA-- 

4* C ALL WATER SPECIFIC HEAT f DENSITY ASSUMED TO BE — 1.C-- EXCEPT HOT WATER 

5* DIMENSION D AP(10) ? XIN(!G),0UT(1G) 

6* DIMENSION X(6) f YC6>,GIN(3) f XNTU(6>fEFFN(61»PFR(&) 

7# C 

B* C — METR IC~L*BRITISH UNITS USED. JWRITF=1 WRITE ALL,JWRITE=0 NO WRITE 

9* C KLBHR-otGPM FOR FLOW INPUT. KLBHR=l f LBS/HR INPUT 

1C* DATA METRIC/C/, KLBHR/l/ f JWRTTE/l/ 

11* DATA METRIC/C/,KLBHR/ 1/,JWRITE/Q/ 

12* DATA PDKG/.4S36/ 

13* DATA TFTC1/32./ , TFTC2/1 .8/ 

14* DATA CALBTU/3. 96831/ 

15* C 

16* C C0NST1 l C0NST4 A R F CONCENTRATION LIMITS 

17* DATA C0NST1/O.4/, C0NST4/D.6R/ 

18* C A-B-C-D-CONST ARE LIMITS FOR EVAP., ABSORP., COND.tC GENERATOR 

19* DATA ACON$T/l./ r 8CONST/1.296/ f CCONST/1.423/ f DCONST/1.919/ 

2C* DATA TELO/2.22/, TE2SET/4.43/ 

21* C COP LIMITS — HEAT LOSS FACTOR 

22* DATA COPHI/0.93/ ? C0PL0/C.6G/ f FQG/l •□/ 

23* C EFFNX-D. 71428 FOR T5=135 F EFFNX~(TG-T5 )/( TG-TA1 

24* DATA FFFNX/9. 71428/ 

25* C 

26* C 

27* C 

28* C 

29* I NP U T - 1 

30* 99C CONTINUE 

31* H 1-C . 

32* H5=G. 

33* H7rC. 

34* H 8 - C • 

35* H ID ” 0 . 

36* C 

37* C UA VALUES ARE PER TON BASIS 

38* UAG-456.981 *PDKG 

39* U AC -1011 .869*PDKG 

4 0* UAE rlSC 3. 294* PD KG 

4 1* UAA = 1102 . 8*PDKG 

42* UAX=11S.929*PDKG 

43* UAXO-1 18 .929*PDKG 

44* G SO “14 4 .077*PDKG 

45* G WG - 1 31 .323*PDKG 

46* G HQ -150. 786*500. *0.975 

47* GEQ-417.600 *50D. 

48* 6C 0-553.674*500. 

49* T0NC=174. 

50* C 

5 1* C 

52* C 

53* C 

54* PBR(3)-UAG 

55* PPR { 2 ) - U A C 

56* PER ( 3 ) -U A E 

57* P BR (4 ) - U A A 

5 8* PBR ( 5 ) - U A X 

59* C 

6 0 * JREF-G 

61* 7 CONTINUE 

62* T ON X “ X I N ( 7) /1200D. 

63* T0NREF=T0NQ*0.5 

64* TONM IN= TONO*C : . 1 

65* TONMAy-TONG*! .2 
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6 6 * 

C 




67* 

C 


K L PH R “0» G PM G PM INPUT. — KLBHR=1,BPH 

INPUT FOR FLOW RATE 

68* 



BPH- 500 • 0 


69* 



IF(KLBHR.GT.O) BPH = 1 .0 


7 0* 

c 




71* 


9 

CONTINUE 


72* 

c 




73* 

c 


FLOW RATES ART PER TON BASIS 


74* 

c 


PFR TON RASIS 


75* 



T0N = 1 .0 


76* 



KTONlrO 


77* 



KTON2-0 


78* 



KT0NX--1 


79* 



X 4 - D • 6 7 


8 L * 



X 1 - E • 4 1 


81* 

c 


S --IS A CONTROL COST. — S=0. OUT“INLET — S = 1 OUT CALCULATED 

82* 



S - 1 • c 


63* 


1 1 

CONTINUE 


84* 



GS-GSO 


85* 



KCOP = T 


66* 


16 

CONTINUE 


87* 



IF(COPX.LE.COPLO.OR.COPX.GE.COPHI) 

COPX =C .722 

86* 

c 




89* 

c 




90* 



GH- X I N ( 1 )/T0NX*PDKG*BPH*G.975 


91* 



GC=XIN( 2)/T0NX*PDK6*BPH 


92* 



GfrXINt 3) /TONX*PDKG*BPH 


93* 



TM1= (XINC4 1- TFT Cl) /TFTC2 


94* 



TCl=(XIN(5)-TFTCn/TFTC2 


95* 



TE1-(XIN(61-TFTC1) /TFTC2 


96* 



Of. - 3023.9573 


97* 



TC12-TC1 


98* 



TG-TH1 


99* 



TArTCl 


IOC* 



TC-TC1 


121* 



TE-TE1 


1 u2* 



T3=TA 


10 3* 



T 5“ T6 


1 J 4 * 



DO 6 1=1,7 


105* 



IF( XIM(I ) .LE.0.C001) S = G.G 


106* 


6 

CONTINUE 


107* 



IF(S.LE.u.3D01) GO TO 5 


108* 

c 




iv » * 

c 

ASSUME TUBE SIDE WATER FILM COEF. IS 

CONTROLLING 

11C* 



GH1-(XIN( 1 >/GHb)**D.8 


111* 



GC1“ (XIN(2>/6Cu)**0.8 


112* 



GF1 = tXIN (3 » /GEC )**0.8 


113* 

c 




114* 

c 


TOTAL BASIS IN METRIC UNITS 


115* 



GIN ( 1 > = GH/GH1*T0NX 


116* 



GIN (2)=GC/GC1*T0NX 


117* 



GIN ( 3 ) = GE /GF 1*T0NX 


118* 

c 




119* 



DO 10 1=1,4 


120* 



IF < GINtI ) .LE.C.O I GO TO 8 


121* 



XNTU(J 1=PBR (I 1/GIN (I) *1 • 0* T ON C 


122* 



IF( I .F0.4 1 XNTUCI )=PBR(I)/GIN(?)*TONO 

123* 



IF(XNTU(I>.GE.1C.) GO TO 8 


124* 



EFF N (I ) rl .3-EXPI -XNTU ( I ) ) 


125* 



GO TO 10 


126* 


s 

CONTINUE 


127* 



EFFN(T)=U.9°9 


128* 


1J 

CONTINUE 


129* 

c 




15 0* 

c 




131* 


5 

CONTINUE 


132* 



IF(TIME.LE.TIMED.0R.TE1.LE.TE2SETJ 

TE 1 = T E2SET 
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133* 



I F ( TL l.LE .TE2SET > S=?.D 

134* 



TE2-TF 1-QE/GE*S 

135* 



TH2nTHl-(QE/C0PX)/6H*S 

136* 



TC2-TC1-M 1.3+1.Q/C0PX )*CE/GC*S 

137* 



IF 1 S .LE . O.noOl ) GO TO 80 

138* 

C 



139* 



TF-TE!-<TE1-TE2)/iFFNC3) 

140* 



TG=TH1-(TH1-TH2 >/EFFN ( 1) 

141* 

c 


ASSUMFD A VALUE FOR TA 

142* 



T A2- CTC1 + TC? ) *□ • 5 

143* 



TA-TCI-(TC1-TA2 )/FFFN(4) 

144* 



KTA = E 

145* 


15 

CONTINUE 

146* 

c 


IF ( TA .LF . (TC1 + 1 .0) ) GO TO 43 

147* 



TC=TC2/EFFNf2>-U.Q/EFFN(2>-l.C)*CTCl+EFFN(41*CTA 

148* 

c 



149* 

c 



150* 



I F ( TC.LF.TA ) GO TO 40 

151* 



IF C TF.GE . TA ) GO TO 4 1 

152* 



IF( TC.GE.TG) GO TO 41 

153* 

c 



154* 



XI~(49.04+1 • 125*T A-TF )/( 1 34 .65+3.47*TA ) 

155* 



X4 = <49. 04+1. 125*TG-T 0/(1 34. 650. 47*TG) 

156* 



I F < X ] .LT .C0NST1 ) GO TO 45 

157* 



IF(X4.LE.X1 ) GO TO 43 

158* 

c 



159* 



HP-TC-25 .0 

160* 



H1G-572.8^0.417*TE 

161* 



GR-QE/(H10-H8 ) 

162* 



GS-GR*X4/(X4-X1 ) 

163* 



G W- G S * < X 1 / X 4 > 

164* 

c 



165* 



C X 1 - 1 • D 1 — 1 .?3*X1+0.48*X1**2 

166* 



CX4D1 .01-1 . ?3*X4+0 .48*X4**2 

lb 7* 



G SC 1 - GS *C X 1 

168* 



6WC4 rGW*CX4 

169* 



CPATI0-GWC4/GSC1 

170* 

c 



171* 

c 

ASSUMING original film coef. equal on both sides.- 

172* 



Fl = 2.r 

173* 



f 2- 1 .r 

174* 

c 

FI 

-2 *F 2 - 1 * HG S - HG W • — F lrl,F2=0 f HGS>>HGW.--Fl=2.5,FZ=2 

175* 



RCS - ( (GSv/GS)*( TONO/TONXi ) * * G • 8 

176* 



R GW - ( (GUO/GW )*(TONQ/TONX) ) * *u • 6 

177* 



UAX-F 1 *U A X Q * ( 1 • 0 / ( RGS + F2*RGW ) ) 

178* 

c 



179* 



IF (GWC4 .6T . GSC 1 ) GO TO 17 

180* 



XNTUX-UAX/GWC4* (TONP/TONX > 

181* 



GO TO IP 

182* 

c 



183* 


17 

CONTINUE 

184* 



X N T U X -U AX/GSC1* (TONO/TONX) 

185* 



CRA T IOr&SC 1 /GWC4 

186* 


18 

CONTINUF 

187* 

c 



188* 



IF(ABStl.C-CRATIO) .LT .0.31) GO TO 13 

189* 



IF C ( XNTUX* M .-CPAT IC> ) .GE . 1 T. ) GO TO 12 

190* 



EXPX-EXP( -XNTUX*t 1 .H-CRATI0I1 

19 1* 



E FF N X - ( 1 .O-FXPX )/ ( 1 •j-CRATIO*FXPX ) 

192* 



GO TO 14 

193* 

c 



194* 


12 

CON T I NUF 

195* 



EFFNX-L.999 

196* 



GO TO 14 

197* 


13 

C ON T INUE 

198* 



E FFNX-XNTUX / ( 1 . 0 + X NTUX ) 

199* 


14 

continue 


FOR 


UAX 
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230* 

C 



201* 

C 



2C2* 



T5-TG-EFFNX* (TG-TA > 

2C3* 



T3-TA+(EFFNX*CRATT0*(TG-TA > ) 

204* 

C 



205* 



Hlr (42.81-425.92*X1+404.67*X1**2)+CX1*TA 

2u6 * 



H5r<4?.8l-425.92*X4+404.67*X4**2>*CX4*T5 

2^7* 



H 7-5 72. 8 + 0.4 6 *T 6-0 ,04 3*TC 

238* 

C 



2 u9* 

C 


F06=1.0 f N0 HEAT LOSS. — F QG > 1 . 0 , C < 2 . 0 ♦ X HEAT LOSS 

210* 



IF(FOG.LE .1 .0) F QG- 1 . 0 

211* 



QGr CGW*H5-6$*Hl*GP*H7 )*FQG 

212* 



QC-GR*(H7-H8 1 * ( 1*0 + 06/ (06 + QF )*C 1 • Q“F QG) 1 

213* 



QAr tGW*H5-GS*Hl+GR*Hir >*C1 .0 + CG/(0G*QE >* ( l.D-FOG ) ) 

214* 



COP=QF/OG 

215* 

C 



216* 



TC12=TC1+QA/GC 

217* 



IFCTC12.GE.TC2J TC12-TC1+1.0 

218* 



TAX=TCl-( TC1-TC12) /FFFN(4 ) 

219* 

C 



22C* 


40 

CONTINUE 

221* 



IF ( TC .LE . TA ) T AX-TC 

222* 



IF ( ABS ( T AX-TA 1 .LT . 0. COCCI ) GO TO 41 

223* 



IF IKTA.EQ.53) GO TO 41 

224* 



TA- (TAX+TA ) *0.5 

225* 



K TA-KTA + 1 

226* 



GO TO 15 

227* 


41 

continue 

228* 

C 



229* 



IF ( APS ( COPX -COP ) .L T .0 . :cur 1 J GO TO 42 

230* 



IF (KC0P.EQ.5C) GO T042 

231* 



COPX-(COPX+COP)*3.5 

232* 



KCOP-KCOP+1 

233* 



GC TO 16 

234* 


42 

CONI 1NUE 

235* 

c 



236* 



X (1 )-TE?-TE 

237* 



X (2 ) = TA-TC12 

238* 



X (3 * -TC-TC2 

239* 



X (4 * = TH2-TG 

24 C* 



Y (1 ) - A C 0 N S T 

241* 



Y (2 ) -BC ON S T 

242* 



Y (3 )=CCONST 

243* 



Y <4 > -DC 0 N S T 

244* 

c 



245* 



DO 47 I - 1 ,4 

246* 



I F ( X ( I ) . LT . Y ( T ) . AND .K TONX .E Q . 1 ) GO TO 62 

24 7* 



IF(X(I).LT.Y(I ) } GO TO 45 

248* 


47 

CONTINUE 

249* 

c 



250* 



IF(X1.GT.C0NST1.AND.X4.LT.C0NST4.AND.X4.GT.X1> GO 

251* 

c 



252* 


4 5 

continue 

253* 



IF (KTCN2-1QG) 49,43,43 

254* 


49 

continue 

255* 



tonmin=tonx 

256* 



T0NXr(T0NX+T0NMAX)*C.5 

257* 



ITONX-IFIX(TONX) 

258* 



TONX~FLOAT (IT0NX)+1«u 

259* 



KT0N2-KT0N2+1 

260* 



GO TO 11 

261* 

L 



262* 


43 

continue 

26 3* 



I F ( K T ON 1 - 1 0 C ) 44,5r,E: 

2o4 * 


44 

CONTINUE 

?o 5 * 



TONM AX=TONX 

266* 



T0NX=(T0NX+T0NMIN)*~.5 
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267* I TONX-IFIX ( TONX ) 

26 ft* IFfTON.LE.l.U.AND.TONX.LE.ITONMIN + l.OM 60 70 50 

269* IF( TONX. LE. TON) GO TO 50 

2 7 C* TONX=FLOAT< I TONX) -1.0 

271* KT0N1=KT0N1#1 

272* 60 TO 11 

273* C 

27*4* *46 CONTINUE 

275* IF(KTONX.EQ.l.OR.TONX.LE.O) GO TO 60 

276* IFCTON.GE.TONX) GO TO 50 

277* C CHECK MAX. STRONG SOLUTION PUMP RATE 

278* GSTC-GSO*TONO 

279* GSPUMP=GS*TONX 

280* IF(GSPUMP.GT.GSTO) GO TO *4 3 

281* C 

282* T ON z T ON X 

283* IF(TE2.LT .TE2SET .OR.TE .LE.TELO) GO TO *4 8 

29*4* GO TO *4 9 

285* *48 CONTINUE 

286* TON=TON*( TE 1-TE2SET)/ (TE1-TF2) 

287* C 

288* 50 CONTINUE 

289* IFCKTONX.EQ.l) GO TO 60 

29 D* TONX-TON 

291* IF(TON.LE.l.D) S=O.G 

292* K TO N X - 1 

293* GO TO 11 

29*4* C 

295* 63 CONTINUE 

296* C 

297* C 

298* C ALL SPEC HEAT = 1. 

299* U AG - X NT U ( 1 ) *GH 

3C0* UAC-XNTU<2)*GC 

3C 1 * UAE -XNTUC 3)*GE 

332* U A A - X NT U ( *4 ) *GC 

303* C 

30*4* A - A L OG ( 10.0 > 

315* B-1555.3/(TF+273.15) 

3C6* C - 1 1.2914E4/(TE + 273.15)**2 

307* PE-EXP< A*I7.8553-B-C> ) 

328* B-1555.3/ITC+273. 15) 

339* C~ll .2*4l*4E*4/ITC + 273.15)**? 

31 G* PC“EXP(A*(7. 8553-9 -C) ) 

3 11 * C 

312* 80 CONTINUE 

31 3* QC- OG*C A L B T U 

31*4* QC-OC*CALBTU 

315* OE-OE*CALBTU 

316* OA-OA*CALBTU 

317* H 1-H1*CALBTU 

318* HFrH5*C ALBTU 

319* H7=H7*CALBTU 

32 C * H8-H8*CALBTU 

321* Hlu=HlG*CALBTU 

322* C 

32 3* T ?-T3*TFTC2*TFTCl 

32*4* T5=T5*TFTC2*TFTC1 

325* C 

326* TH1 -TH1*TFTC2+TFTC 1 

32 7* TH2-TH2*TFTC2 + TFTC 1 

328* TC1=TC1*TFTC2*TFTC1 

32 9* TC12-TC12*TFTC2-* > TFTC1 

33 C* T A2-TC12 

331* TC2=TC2*TFTC2*TFTC1 

332* TF1-TE1*TFTC2+TFTC1 

33 3* T r 2-TE2*TFTC2-»-TFTC 1 


I 
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334* 
335* 
336* 
337* 
338* 
339* 
340* 
3*1 1 * 
342* 
343* 
34 4 * 
345* 
346* 

34 7* 
348* 
349* 

35 Q* 
351* 
352* 
353* 
354* 
355* 
356* 
357* 
358* 
359* 

36 G* 
361* 
362* 
363* 
364* 
365* 
366* 
367* 
368* 
369* 

37 G* 
371* 
372* 
373* 
374* 
375* 
376* 
377* 
378* 
379* 

38 D* 
381* 
382* 
38 3* 
384* 
385* 
386* 
38 7* 
388* 
389* 
390* 
391* 
392* 
393* 
394* 
395* 
396* 
397* 
398* 
399* 
400* 


TL -IE *TFTC2+TFTC1 
TA =TA *TFTC2+TFTC1 
TC -1C *TFTC2+TFTC1 
TG rTG *7F7C2+TFTC1 
C 

UAG-UAG/PDKG 
UAC-UAC/PDKG 
UAE-UAE/PDKG 
UAArUAA/PDKG 
UAXrUAX/PDKG 
6H-GH/PDKG 
G C- GC/PDK G 
G A ” GC 

GE-GE/PDKG 
GP=GR/PDKG*S 
6V-GW/PDKG*S 
GS“GS/PDKG*S 
GSC 1-GSC1/PDKG*S 
GWC4-GWC4/PDKG*S 
C 

QGT -QG*T 0N*S 

OCT “ OC*T 0 N* S 

QETrTON 

OAT =OA*TON*S 

UAGT=UAG*TON 

U AC T -UA C * T 0 N 

UAET=UAE*TON 

U A A T =IJA A * T 0 N 

GHTrGH*T0N/8PH/0.975 

GCT -GC*T0N/BPH 

gat=gct 

GETrCE*TON/PPH 

DT1?E=TF1-TF2 

DT12A-TA2-TC1 

DT12C-TC2-TA2 

DT12G=TH1-TH2 

OTE2=TE2-TE 

DTA?=TA-TA2 

DTC2=TC-TC2 

DTG2-TH2-TG 

C 

C 

IF( JWRITE .EO.O) GO TO 59 

IFI AMODf TIME, 1.000) •GT.DELT) GO TO 500 

WPITE<6,95) 

T ON XIN-XIN(7> 

X IN ( 7> = T0NXIN/12Q0G* 

WRITE (6,98> KTA ,KC0P,KT0N? ,KT0N1 , CXlN (I > »Irl t 6),T0NX,XIN<7),lNPUT 
X IN t 7 ) - T 0 NX I N 
WRITE (6,403) 

WRITE (6,4 02 I XI ,X4 ,CX 1 ,CX 4 ,GR,GS,GW,GSC1 , GWC 4 , E X PX , CRA T 1 0 , U A X 
WRITE (6,407) 

WRITE (6,402 ) TA f T5,T3,TG,Hl,H5 f H7,H8 f H10,XNTUX,EFFNX,COP 
WRITE (6,406 ) 

WRITE (6,4 05 ) GE , TE 1 ,TE?,TE ,XNTU(3>,EFFN<3 ) , P E , li A E , Q 5 , G E T , UAE 
XT , OFT 

WRITE (6»4Q4 ) G A , T C 1 , TA 2 , T A , XN TU ( 4 ) ,EFF N (4 > , P E , U A A # 0 A , GA T , U AA 

XT , OAT 

WRITE (6,401 ) GC,TA2,TC2,TC ,XNTU(2),EFFNC2),PC,liAC,0C,GCT,UACT,0 
XCT 

WPITE(6,400) GK , TH 1 ,TH2,TG ,XNTU(1),FFFN(1),PC,UAG,0S,GHT,UAG 

X T ,QGT 

400 F0RMAT(1X,4HG ,7F1C.3,5E10.3//) 

4 01 F OR MAT(1X,4HC , 7F 1 Q . 3 , 5E 1 D • 3/ ) 

402 FORMAT( IX ,4HX ,!2FlC.3/> 

403 FORMAT(5X,120H XI X4 CXI C X 4 GR 

X GS GW GSC 1 GWC4 EXPX CRATIO UAX 
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401* X ) 

MG 2* 4 04 FORMAT ( IX ,4 HA , 7F 1C . 3 ,5E10 . 3/ > 

4 J 3 * 4 05 FOR MAT (IX , 4 HE *7F 1G .3 f 5E 10. 3/ ) 


4 04 * 

406 F0RMATC5X r 120H 

G 

T 1 

T 2 

T 

NTU 

405* 
4 06* 

X EFFN P 

X > 

UA 

0 

GT 

U A T 

OT 

407* 

407 F0RMATC5X ,120H 

T A 

T 5 

T 3 

TG 

HI 

4 l. 8 * 

X H 5 H7 

H8 

H10 

NT UX 

EFFNX 

COP 


4j9* X ) 

410* C 

411* WRITE (6»961 

412* WRITE (6,97) DT12E,DTE2,TE2 t DT12A f DTA2,TA2,DTl?C,DTC2,TC2 f QT12G f DTG 


413* 


X 2 * T H2 



414* 


95 FORMAT (5X * 1 2 7 H K T A KCOP 

K TON 2 KT0N1 5HT 

-GPM 

415* 


X GCT-GPM GET-GPM TH1-F TC1-F 

TE1-F TON-CAL 

TON-ST 

416* 


X A RT NO. ) 



417* 


96 FORMAT CSX , 1 2DH DT12F DTE2 

TE2 DT12A 

D T A 2 

418* 


X T A 2 DT12C DTC2 TC2 

DT12G DT 62 

TH2 

419* 


X ) 



420* 


97 F OR M AT ( IX , 4 HTEMP ,1 2F10.3/// ) 



421* 


98 F0RMAT(5X,4 (18,2X1 ,8F10.3*3X f I3/) 



422* 


99 F0RMATC7F1Q.1 J 



423* 

C 




424* 


59 CONTINUE 



425* 


5 HO CONTINUE 



426* 


OUT ( 1 ) “ T H 2 



427* 


OUT (2>=XIN( 1) 



428* 


OUT ( 3 )=TON*OE 



429* 


OUT (4 > = TH2 



4 3 C* 


OUT (5 >=OGT 



431* 


OUT (6>=TF2 



432* 


OUT (7)~XIN(3) 



433* 


OUT (8 )=TC2 



434* 


OUT (9)-XIN(3) 



435* 


OUT ( 10)rTON*OE 



436* 

C 




437* 


INPUT- 1+INPUT 



438* 


R F T URN 



439* 

c 

STOP 



44 D* 


END 
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Figure 1. - Flow diagram of single-stage LiBr-HgO absorption unit 
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